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METHOD OF DESIGNING CASCADE BLADES WITH PRESCRIBED VELOCITY 
DISTRIBUTIONS IN COMPRESSIBLE POTENTIAL FLOWS 

By George R. Costello 


SUMMARY 

By use of the assumption that the pressure-volume relation is 
linear, a solution to the problem of designing a cascade for a 
given turning and with a prescribed velocity distribution along 
the blade in a potential flow of a compressible perfect fluid was 
obtained by a method of correspondence between potential flows 
of compressible and incompressible fluids. The designing of an 
isolated airfoil with a prescribed velocity distribution along the 
airfoil is considered as a special case of the cascade. 

If the prescribed velocity distribution is not theoretically 
attainable, the method provides a means of modifying the dis- 
tribution so as to obtain a physically significant blade shape. 
Numerical examples are included. 

INTRODUCTION 

In order to control boundary-layer growth, transition, and 
separation in the design of a cascade for a given tu rning , it is 
advantageous to prescribe the velocity along the blade as a 
function of the arc length along the blade and then to com- 
pute the blade shape. For an incompressible fluid, several 
solutions to this problem have been obtained. (See refer- 
ences 1 to 3.) 

A similar solution for the two-dimensional potential flow of 
a compressible perfect fluid was developed at the NACA 
Lewis laboratory during 194S-49. This solution is based 
on the assumption that the pressure-volume relation is 
given by a linear approximation to the isentropic curve 
instead of the true curve. The flow pattern of the com- 
pressible fluid is obtained by a transformation from a cor- 
responding flow of an incompressible fluid using the trans- 
formation developed by Lin (reference 4). 

The method of solution consists in using the free-stream 
velocities upstream and downstream of the cascade and the 
prescribed dimensionless velocity distribution along the 
blade to select a suitable incompressible potential flow about 
the unit circle and then to determine the mapping function 
that transforms this incompressible flow into a compressible 
flow about a cascade of airfoils. The image of the unit 
circle under this mapping gives the cascade with the pre- 
scribed velocity distribution along the blade, provided the 


velocity distribution is theoretically attainable. If the 
velocity distribution is unattainable, methods are given for 
modifying the distribution so that a physically significant 
profile is obtained. 

The problem of designing an isolated airfoil with a pre- 
scribed velocity distribution is considered as a special case 
of the cascade. 

THEORY OF METHOD 

CASCADE 

In reference 4, Lin shows that if the pressure-density rela- 
tion is (Symbols used in this report are defined in the 
appendix.) 

P=C,-^ ( 1 ) 

then the compressible potential flow about a cascade of 
blades can be obtained by transforming the incompressible 
flow about the unit circle in the following manner: The com- 
plex potential function FIX) for the incompressible flow due 
to two complex sources at f=a l and f—a> outside the unit 
circle |f|=l is 

F(g)=A log, (r —af}+A log. ($--=-) -KB log, (f— 

3n °s*0-r> z ' < 2 > 

where A and B are complex constants with Be A> 0 and 
Re A=—Re B, and D is an arbitrary complex constant. 
The bar indicates the complex conjugate. The mapping 
between the s-plane and the Kplane defined by 

dz = £f(f)(r-ai) -I (r-a 2 )- 1 [F v G‘)l 2 [s)f)l _1 (f ~ a-iXf- a 4 l df 

(3) 

gives a compressible flow with the linear pressure-volume 
relation past a straight cascade of identical blades in the 
s-plane with the velocity potential $ c and stream function 
given by 

4‘c J ri'Pc=F(g) (4) 
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provided that g(j) is chosen to satisfy the following re- 
quirements : 


(a) The function p(f) is regular in closed region 2 ?' 
defined by |f| > 1. 

(b) The function g(£) ^ 0 in 22 , except possibly at 

one point on the circle where 2 fI, (f)=0 (the order of 
the zero not to exceed 1). ► 

(c) Along the circle |f[=l ^ dz = 0 . 

(d) The function </(f) satisfies the inequality 

f[^(f)]to(f)]- 1 (r-a0(f-o^|<2 in R. J 


( 5 ) 


The magnitude q and the direction a of the dimensionless 
velocity at any point in the 2-plane are given by 


from the given conditions. The constants A and B are 
obtained from the upstream and downstream velocities and 
the circulation and then a y and a 2 are determined by the 
range of the potential on the airfoil. 

Circulation and cascade spacing. — The magnitude of the 
prescribed dimensionless velocity g along the airfoil is given 
as a function of the arc length [g=g(s)] where the total arc 
length is taken to be 2 t and is measured from the trailing 
edge along the lower surface. If Q(s) is defined by 


Q(s) =—q(s) 0 <s<s tt 

Q( s )=q(s) s„<s<2tt 


( 8 ) 


where $„ is the leading-edge stagnation point, then 


2g g - ia ^'(j-Xr-giXr-fflz) 

l + Vl+2 2 0(f) 


( 6 ) 


In order to use this transformation in designing a blade 
with a prescribed dimensionless velocity distribution along 
the blade in a cascade, the prescribed conditions -are used to 
select a suitable incompressible flow about the unit circle 
and to determine the function p(f). 

The prescribed conditions are the velocity distribution on 
the airfoil, the upstream velocity q\i ai and the downstream 
velocity g 2 e f “ 2 . The upstream and downstream velocities are 
related by the isentropic-flow equations with 7= — 1. This 
relation is 

g 2 a cos 2 a 2 __ gi 2 cos 2 ^ , . 

l + g 2 2 ~ 1 + g! 2 ^ ' 


where the axis of the cascade has been taken along the y-axis 
for convenience and the flow is from left to right. (See fig. 1 .) 





Fioukb 1.— Cascade in 2 -plane. 



FLOW IN CIRCLE PLANE 


The flow of an incompressible fluid about the unit circle is 
selected by determining the constants A, B, a h and « 2 in the 
complex potential 


F(t) = A log, (f—aJ+A log, (V -=-) -f- 
B log, (f-a 2 )+5 log, 


Q(s) 


ds 


( 9 ) 



( 10 ) 


The circulation and the spacing of the cascade are related 

by 


<2 


r, 

2i sin cci— g 2 sin a 2 


( 11 ) 


where d is the spacing. The quantities r„ q u g 2 , a u and 
a 2 are known so that the spacing is determined. 

Determination of A . — The value of d from equation (II) 
is used to evaluate A and B because the spacing is also 

given by the absolute value of ^ dz taken along a path 

around a x or a 2 . (See fig. 2.) The axis of the cascade has 
been taken along the y-axis so that 



Figure 2.— Paths oC Integration in {"-plane. 


The second equality comes from the fact that the residues 
at infinity of 


( 2 ) 


0(r)(?— «i) -1 (f— a») _1 
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and 


[F r m 2 ig (fli-'cr-aocr-aj 


in the expression for dz in equation (3) are zero and 
consequently 

■(f) dz=0 (13) 

Jc J ®1 JO 2 

where c is the unit circle. But by equation (5c), 

§ c dz= 0 


so that 


dz 


(14) 


The evaluation of equation (12) in terms of the potential 
F(£) is 


■id=<D 77 — ^ [^COUCf °iKL a A dt 
(r — aO(r — aa) 4j*i ff(r) 


But 

(r) f-a t + l_ 

* ay a 2 

so that equation (15) reduces to 


(15) 

(16) 


< 1S > 


At f =a u equation (6) becomes 


2gi „-fa t _ A(ay — (h) 

1 + Vl+Si 2 ff(«0 


(19) 


which on writing 


reduces to 


K r 




l + Vl+2i 2 
g(a J _A 


(ffi — a 2 ) A 1 


e ! “i 


( 20 ) 


Substitution of the values from equation (20) in equa- 
tion (18) gives . 

id =27ri [(^+^)« Iai ] (2 1) 

Hence, the bracketed expression in equation (21) must be 
a real number and 


4 A+AKy- 

— 4 Ky ^ 


(22) 


where 




From equation (22), 


or 


(4 — Ky)ImA 

(4 +KfiBeA tan ai 

A.-Z-K- 

ImA= — 1 * i?e A tan a t 

4— IsTr 

ImA — — V 1 + 2 ! 2 -He A tan a t 


(24) 

(25) 


which gives the relation between Be A and In A. 

Substitution of the value of Im A from equation (24) in 
equation (23) yields 

Hwy Be ’ A+ (wy s ‘ ,a *“■ “■ 

/4iZA 2 n,, , 

=( — — j Re 2 A sec 2 ar 

’ 4 +KC- 


or 


4JEi 


-i?e A [sec ad 


Hence 


d=2w(^F^ReA [sec ad) (26) 

Substitution of the value of d from equation (26) in 
equation (11) gives 

2x ( iKy ReA ^ sec ai 0 = ^r 


or 


ReA-- 


sm ay — q 2 sm a -2 
'kKyTc [cos a t | 


27r(4+AT 1 2 )(g l sin ay—q_ 2 sin a^ 


(27) 


and Re A is now determined. By use of this value of Re A 
in equation (25), Im A is obtained. Hence A is completely 
given by equations (27) and (25). 

Determination of B . — From equation (12) 
id=— dz 

= ~^ 2 ff(f)(f— ®i) -I (r— aJ _I <ff- + 


j <j^ dtr 


id=— 2tti 


• 9(fl2) , 1 B 2 (g2— ad 27ri 


(a 2 — a0 4 


9(0-2) 


„ . Be ia 2 2-ri , 
= — 2iri — BK 2 e ,a o 

xv 2 4 

_ . r/4B+i£, 2 5\ . 1 

=- 2,n L(-^— > 2 J 


where 


K 2 = 


2g 2 


1+ Vl+2? 
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The bracketed expression must be real, so that 


(4 —K 2 2 } Im B 
(4+iTs 2 ) ReB 


tan 0:2 


But 


ReB= —Re A 


and equation (28) can be written 


(28) 

(29) 


4>(8) = — 2ReA tanh -1 ^ -\-(Im A+Im B) 


. tan 8 
t£m fTiWT 


(Im A - 


- Tm. 


Im B—^^j^Re A tan a% 

— Vl + ff 2 2 Re A tan Qto (30) 

Consequently, B is determined by equations (29) and (30) 
because Re A is known from equation (27). 

Determination of a t and a 2 . — After A and B are known, 
the po'ints a t and a 2 are to be selected to satisfy tlie single 
condition that the range of potential on the circle must equal 
the range of potential on the airfoil, that is, 

M 2 tt) - 4> c (s n )=F[e i <?t + 2 -)] -F(e ".) (3 1) 

where 8 t and 8 a are the trailing-edgo and leading-edge stag- 
nation angles, respectively. This condition is only one 
equation in two unknowns, a, and a 2 ; consequently, the values 
of a l and a 2 are not uniquely determined. By imposing an 
additional restriction that a x and a 2 are real and . 


ai =—a 2 (32) 

unique values are obtained for a x and a 2 in all cases. In par- 
ticular problems, however, some other restriction may be 
more useful, such as assigning a definite value for a x and 
computing a 2 . 

With the restriction given in equation (32), it is possible 
to express 8 t and 8 n in terms of a x and substitute these values 
in equation (31), but the resulting equation cannot be solved 
explicitly for a } . 

One method for obtaining a x is as follows: Let 


Oi=e* I 


a 2 =—e k 


j 


(33) 


where 7c > 0 . Then equation (2) becomes 

F(X)=A log e (r— e*)+ J log e (r-e“*)+-8 log«(rH-e*)4- 
B log, (f ■ +e~*)+D 

or, for points on the unit circle f =e (e 

i°g. 

ilm A log, I Urn B log, (34) 


where D has been chosen to make <jn(8 t ) = 0 and the angle 
convention is 

t , . _. sin 8 

-2 iIHCT<2 

and tan -1 is taken in the same quadrant and same 

direction as 8. 

The velocity on the circle v(8 ) is 


>( 6 ) 


2 Re A sin 8 cosh k ( ImA+ImB) tanh k sec 2 0 
~ cosh 2 k— cos 2 8 ' tank 2 /c+tan 2 6 

,T a r t>\ sink k cos 8 

(ZmA - ImS) S3n+S?e - • ■ 


cosh 2k — cos 28 


[2 Re A sin 0 cosh k+ 


(Im A — Im B) cos 8 sinh k + (Im A+Im B ) sinh k cosh k] 

(30) 

Equation (36) can be further simplified by defining X as 


. (Im A—ImB ) sinh k 

tan X= 2 IferiL coshT 


(37) 




Then 


, . 2 V 4: Re 2 A cosh 2 k 4 -(Im A—Im B) 2 sinh 2 k 

~~ cosh 2k— cos 28 




sin 8 cos X + 


„ . , . (ImA-\-ImB) sinh k cosh k ~\ 
cos 6 sm X 4- ■■■ ,„. - ■ , ■■■ == — • ? -• « 

■yj4:Re 2 A cosh 2 k+(Im A—ImBf sinlr/r J 

4:ReA cosh k sec X F . , . , . (Im A+Im B) sinh /l*~1 

cosh 2t-c ' os 2 8 L 8m (S+X)+ ®T5S J 

(38) 


The stagnation angles 8 t and 6 n are therefore the roots of the 
equation 

—(Im A+ImB) sinh k 
2 Re A sec X 


sin (0+X) = 


(39) 


The desired value of k is obtained as follows : 

(1) Assume a value of k. 

(2) Compute X by equation (37). 

(3) Obtain 8 t and 8 n from equation (39). 

(4) Compute <t> i (8 l +2'ir) — cb t (d n ). 

(5) Repeat (1) to (4) several times to obtain a plot of 
0i(0<+2x-)— <j>i(d„) as a function of k. 


Hence, <£(0) may be written in the form 
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( 6 ) Interpolate to obtain k such, that 

4 h(Pi "r 2ir) — 4>t{8n) — — $>e(0 

With k determined., the flow about the circle is known. 
The potential 4 h{ 8 ) and the velocity v(6) for points on the 
circle are given by equations (35) and (38), respectively. 

function g(f) 

The function git) can be computed for points on the unit 
circle by using the prescribed velocity on the airfoil and the 
velocity on the unit circle to determine the real part of git). 
The imaginar y part of 17 (f) can then be computed by Poisson’s 
integral. Because of the restrictions imposed by the given 
conditions, however, g(t) is actually obtained in a slightly 
different manner, as shown in the following sections. 

Airfoil with pointed trailing edge. — If an airfoil with a 
pointed trailing edge is desired, then git) must vanish at the 
trailing-edge stagnation point f— Hence, git) can be 
written in the form 


ff(r)=(l-y)V«» 


(40) 


where/(f) is regular in the exterior of the unit circle and 


— !-i 

X 


(41) 


where 5 is the included trailing-edge angle of the airfoil. 
(See reference 4.) ‘ 

Values of g(±e*). — Because the velocities are given for the 
compressible flow upstream and downstream of the cascade, 
the value of git) at f=±fi 6 is determined from equation ( 6 ), 




gi— <?) 


2 

K, 




— 2 e k Be ia * j 

: k 2 J 


(42) 


In order that git) have these values, /(f) is written in the 
form 


where 


fit)=Cit)+ [e2k n p ~ 2 - e ~ - } Bit) 

g(D -4 (1 +1) loe. [g ^$ giy]+ 

1 A e ks \-> r —2 Be. ia *+ k *1 

2 v 1 fJ^L^a+c^-^J 


(43) 


(44) 


and Hit) is regular in the exterior of the unit circle with 


lim fff(f )=0 

— yco 


(45) 


l -r)‘ 


cay 


I s 


H a) 


(46) 


and g(t) will be known when H(t) is determined. For the 
actual computation of the blade shape, only the values of 
git) on the unit circle are needed. Hence, it is only neces- 
sary to compute Bit) for points on the circle. If desired, 
the values of Hit) for any point in the exterior of the circle 
can be obtained from the values on the circle by Poisson’s 
integral. 

Determination of Re H on circle. — By equation (4), the 
potentials <f>ds) and <f>i(6) are equal at corresponding points. 
Thus, by matching these potentials a correspondence is 
established between points along the airfoil arc and the 
circle angles; that- is, s=s{9). By use of this correspondence, 
the magnitude of the prescribed velocity along the airfoil 
is obtained as a function of the circle angle g =g( 0 ). Hence, 
by ta king absolute values of equation ( 6 ) , 

2 g j6) |F’(e i 8 )(e i 6 -^)(e i 9 -ke*)| 

l+Vl+g(fi ) 2 \gie ie ) | ^ 

for points on the circle. Substitution of the value of g(f) 
from equation (46) with f — e te and replacing F'ie is ) by the 
velocity vie) (equation (38)) on the circle give 

2 g( 0 ) 

i+Vi+<z( 0) 2 

e fc jp(g )[(2 cosh 2k— 2 cos 2d) 2 

[2 — 2 COS (g ( _g)]V e t« e C r (e I ' s ) + C2cosh2t-2cos2fl)neHC<!'' <1 )] 


or, with the equation solved for Re Hie iS ) , 


log* 


ReHie iS )=- 

where 


[p(g)|(2 cosh 2k— 2 cos 20/ 
Kje)[ 2—2 cos ie t — e)]2 


>-ReCie tB )+k 


m 


2 cosh 2k— 2 cos 26 


2 g(fl) 

i+VT+W 


(48) 


(49) 


Restrictions on Re Hie iff ). — Equation (45) imposes restric- 
tions on the values of Re Hie iS ), as shown by writing Hit) 
in the form 

hi , h 2 . 


Hit) — AoTyvpT - - - 

where h 0 , hi, h 2 , . . . are complex constants. 

For points on the circle, equation (50) becomes 

Hie ilF )=Re Hie ie )+i Im Hie™) 


(50) 


The restriction on Hit) imposed by equation (45) is neces- 
sary so that /(f) (equation (43)) will be regular. By use of 
equation (43), git) is expressed as 


=Re hnd-'ytRe h } cos jd+Im h s sin jd)+ 

i [Im iin+yi/m h } cos j6—Re h } sin j6)] (5 1) 
j'=i 
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Equation (51) is a Fourier expansion and 


Re A 0 =J- P T Re H(e iS ) d8 

47t Jo 

(52) 

Re Ai=i f Re H(fi u ) cos O dd 
K Jo 

(53) 

Im hi= - f Re H(e te ) sin 0 dd 
* Jo 

(54) 

But equation (45) requires that 


Re h 0 =Im h 0 —Re hi=Im h{=0 

(55) 

Consequently, Re H(e w ) must satisfy the equations 


7 '2ir 

1 Re H(e te ) d8= 0 

(56) 

Re H{e u ) cos 8 dd— 0 

(57) 

I*2x 

I Re H(e i0 ) sin 8 dd=0 

(58) 


Adjustment of Re H(e iB ). — If tlie values of Re H(e u ) from 
equation (48) do not satisfy equations (56), (57), and (58), 
the values must be adjusted until the conditions are 
satisfied. One method for adjusting the function is to 
define Re H{e n ) by 

Re &(e il )=Re H{e iB )~-Re ho— Re hi cos 8—lm hi sin 6 (59) 

where Re h 0 , Re h u and Im hi are given by equations (52), 
(53) , and (54) , respectively. The modified function Re S (e te ) 
will then satisfy equations (56), (57), and (58). This method 
of modifying Re H(e ie ), however, changes the velocity distri- 
bution all along the profile and, if the correction terms in 
equation (59) are not small, these changes in the velocity may 
be extensive because 

2 ? - - 

l + Vl+f 

2g ^ (JRe Ao+Re A 1 cos S+Im Aj sin 9) (2 cosh 2k — 2eos 21) 

i+VhF? 

where 1 denotes the new Velocity. In some cases, conse- 
quently, Re H(e lt ) can best be adjusted to satisfy the require- 
ments by adding to Re H(e ie ) odd and even functions that 
have nonzero values only in small neighborhoods of the points 
0=0 and 0=— 7 r. The particular functions to be added to 
Re H(e tB ) and their range of values depend on the specific 
problem; no general method can be given for determining the 
functions. 


Determination of Im H(e is ). — After Re H(e 10 ) satisfying 
equations (56), (57), and (58) is obtained, the function 
Im H(e K ) is given by Poisson’s integral (reference 5) 

Im H(e ie )=~ f* Re H{e ir ) cot (~^) dr (60) 

where the constant term in Poisson’s integral has been taken 
as zero so that 

Im h 0 =~- f ImH(e ts ) d &— 0 
*rr Jo 

as required by equation (55). Hence, II(e ie ) is determined 
for points on the unit circle by 

H(e tv )=Re H(e is )+i ImII(e iB ) 

Adjustment of g(e id ). — By use of these values of II(e t(r ) in 
equation (46), p(f) is determined for points on the circle, 

g(e ier )= [1 — e { -®>] “ e tcO")+(2 cosh 2A--2 cos 2 0) /*(«< 9] (g 1 ) 

Because of the adjustments in H(e i0 ), g{e u ) may no longer 
satisfy condition (5d). If g(e ie ) does not satisfy the inequal- 
ity for points on the circle, then the values of 17(f) can be 
adjusted to satisfy the inequality by changing the second or 
higher harmonic terms in H(f) or by other methods. If the 
prescribed conditions are theoretically attainable, however, 
then no modification is necessary, not even in Re II(e ie ). 

BLADE COORDINATES 

By use of the values of g{e iB ) that- satisfy all conditions, the 
blade coordinates are obtained from integration of equa- 
tion (3), that is, 

J gm-eY 1 ^ J ^ - (f = ^ (r - +e * ) ^f 

-i( 

which, on replacing F'(e ier ) by v(8)e \ 2 ) and writing 

g(e u ) — gi(0)e io 2i<>) 

reduces to 

3 = J[p 1 ( 0 )(e--^)- 1 -^| e *[*+ 7 +in »] dO 

(62) 

COMPUTATIONAL PROCEDURE FOR CASCADE 

An outline of the procedure for computing the blade shape 
is as follows : 

(1) Obtain ) and T c from equations (9) and (10), 

respectively. t 

(2) Compute Re A, Im A, and Im B by equations (27), 
(25), and (30), respectively. 

(3) Obtain k as outlined in the text. Compute 4n(0) and 
v(6) by equations (35) and (38), respectively. 

(4) Plot 0 e (s) and 4h(6). By comparing the abscissas for 
equal values of these potentials, obtain s as a function of 0, 
which permits writing the prescribed velocity g as a function 
of 8, 2=2(0)- 
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(5) Compute Re H(e u ) by equation (48) and determine 
Re h 0 , Rek u and Im h t by equations (52), (53), and (54), 
respectively. If these values are not zero, then adjust 
Re H(e u ) either by equation (59) or by addition of functions 
so that Re H(e w ) satisfies equations (56), (57), and (58). 

(6) Obtain Im H(e is ) by equation (60), using the adjusted 
values of Re H(e ie ). 

(7) Obtain g(e ie ) by equation (61). The function g(e ie ) 
must satisfy inequality (5d) for points on the circle. If 
g(e w ) does not satisfy the inequality, adjust g(e K ) as sug- 
gested. 

(8) After g(e ie ) has been adjusted to satisfy all conditions, 
the blade shape is obtained by integrating equation (62). 

ISOLATED AIRFOIL 

The problem of designing an isolated airfoil with a pre- 
scribed velocity distribution along the airfoil in a com- 
pressible potential flow with a given free-stream velocity can 
be considered as a degenerate case of the foregoing cascade 
problem in which the upstream and downstream velocities 
are equal and the spacing of the cascade becomes infinite. 
In this case, the singular points f=a L and V~a 2 move to 
infinity and the complex potential F(f) for the incompressible 
flow about the unit circle (equation (2)) becomes 

-F(r)=— ^f+^+^^log t+D (63) 

where V is the incompressible free-stream velocity and is in 
the direction of the negative 5-axis. The mapping function 
(equation (3)) becomes (reference 6) 

dz=g (r) dr-^[F'(r)i%(f)l- l dr (64) 

where g(f) satisfies the following requirements: 

(a) The function g(f) is regular in the closed region' 

R defined by [f |>1. 

(b) The function £/(f)#0 in R, except possibly at 

one point on the circle where F'(f )= 0 (the order of 
the zero not to exceed 1). „ 

(c) Along the circle |f|=l, dz=0. 

(d) The function g(f) satisfies the inequality 

|[^’ , (f)] [ff(r)l -1 l<2 in R 

The velocity potential <p c and the stream function ip e in 
the z-plane are given by 


In order to use the preceding transformation, the prescribed 
velocity distribution along the airfoil g=g(s) and the free- 
stream velocity g,e !ai are used to select the flow about the 
unit circle and to determine the function 5(f). For the actual 
computation of the airfoil, only the values of the functions 
for points on the unit circle are needed. 

PLOW IN CIRCLE PLANE 

For points on the circle the complex potential 

(equation (63)) reduces to 

4 >i(8)=- 2F(cos S-cos 8 t ) (68) 

when D is so chosen that <£ t -=0 at the trailing-edge stagnation 
point t=e ,e ‘ (fig. 3). The velocity is 

v(8)— — ie~ is ^2V sin (69) 

Both 4,1(0) and v(9) are completely determined when r i; V, 
and 6 t are known. These quantities are obtained from the 



prescribed velocity distribution by determining 4> c (s a ) and T c 
from equations (9) and (10), respectively. By equation (66), 
the potentials are equal at corresponding points so that 

(pi, min = min — Vc(^' r.) ) 

[ (70) 

r,=r e ) 

Then 8, is given by (reference 2) 


<k+*> e =F(r) (66) 

The magnitude g and the direction a of the dimensionless 
velocity in the s -plane are given by 


cot 8 t + 8 t = — 


ir 

2 " 




Ti 


T T<t> c (S„) 

2 r c 


and V is given by 


(71) 


H 

1+ Vi+ g 2 


e 


„ im 

ff( f) 


(67) 


V 


Vi 

4v sin 8 t 


Vc 

4 v sin 8 t 


(72) 
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and then 


The flow about the circle is obtained by using these values of 
Ti, V, and 6 t in equations (68) and (69). 

FUNCTION 0(f) 

As in the cascade case, if an airfoil with a pointed tail is 
desired, g(£) must vanish at the trailing-edge stagnation 
point . f =e i9 ‘. Hence, g(f) can be written in the form 




(73) 


where /(f) is regular in the exterior of the unit circle and 


71=1-- 

TT 


(41) 


By using the prescribed velocity distribution on the airfoil 
and the circle velocity, the real part of /(f) can be computed 
for points on the circle and the imaginary part of /(f) can 
then be obtained by Poisson’s integral. Hence g(f) would 
be determined. Because of the restrictions on /(f), this 
computational procedure is modified as shown in the next 
sections. 

Kestrietions on /(f). — Because the free-stream velocity is 
given for the compressible flow, the value of gr(f) at infinity 
is determined by equation (67) 


lim ff(f)= 




K x 


(74) 


where 


By writing 


where b Q) by, b 3 , 


and 


Ky 


_ 2gi 


i+ Vi+<zi 2 

/(?)== &o+4r+'r|+ ■ • 


. are complex constants, then 
lim g( f)=e 4 ° 

r — >oo 

V' 

Re (6 0 )=loge^ 

Im (b 0 )=ai+ir, 


(75) 


Hence, the desired free-stream velocity will be obtained if 
b 0 satisfies equation (75). 

Further restrictions are imposed on /(f) by the requirement 

<J> rtt) <*r-j <j> F'(r)]%(f)]- 1 df=o (76) 

When the residues of these functions at infinity are con- 
sidered, this requirement can be expressed in a more useful 
form; The function g( f) may be written as 


g(X)—e 6 » 


, , by— ne i0 ‘ . . 

1 i — h — z r terms 


in jj> j> 2 J 


^ £f(f) df =2rie b o(jb l — ne i5 ‘) (77) 


From the incompressible flow about the circle, 

2 Vi sin 0 t . F 


F'( f) F- 


and 


[^(fM'toO-)]- 1 -^ (v 2 - 


4V 2 i sin O t — nV a e u ‘+byV* 


terms in 


F» 2 ) 


therefore 


[F'(f)] 2 [g(f)] -1 df =— 2vie- 6 o(4F 2 i sin 0 t —nV*e i8 ‘+biV*) 

(78) 

Then, by use of equations (77) and (78), equation (76) 
becomes 


2irie 6 o (by— ne i6 ‘)-\-- 2Trie~ b oV 2 (4isin 9 t —ne i0 ‘+by)—O 


or 


V 2 _ 

2irie~ b ° [e 2Reb a (b l —ne ,0 ‘) — ^-(— 4i sin 8 l —ne~ il> ‘-\-Oi))=0 

(79) 

Equating the real parts of equation (7 9) gives 
^e 2Reb °-^-^(R e by—ncos 6 t )—0 
V 2 

But e 2Reb o — cannot be zero except for zero free-stream 
velocity, hence 

Re by=n cos 8 t (80) 

The imaginary parts of equation (79) give 

V 2 

e 2 Reb 0 ^— 72 , sin 0 f )-f— — (Jm by— n sin 8 t )~ f-F 2 sin ff t =0 


or 


Im by=n sin 6 t - 


F 2 sin 9 t 

T/2 

e 2 n«6 0 u.J_ 
4 


But from equation (75) 


F 2 

«2 Rtb 0 — _ 

^ ~Ky 2 


hence 


Im by= sin 6, 


(81) 
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Consequently, equation (76) will be satisfied, if Re &i and Ira 6 t 
satisfy equations (80) and (81), respectively. 

For convenience in computation, these restrictions on the 
values b n and b l are transformed into restrictions on the 
values of Re j(e u ) . For points on the circle, 


f(e ilr ) =Ref(e iB ) + ilm f(e ls ) = b 0 +~+^- e - 


=Reb 0 -\-^P,(Re b 3 cos jd-\- Ira b 3 sin j8)+ 

/=i 


Equation^(82) is a Fourier expansion and 

r 1 r*‘ +2r 


'—Re bj sin. j0)J 

(82) 

and 


\e is )dd 

(83) 

) cos 8 dd 

(84) 

) sin 8d8 

(85) 


Substitution of these values of Re b 0 , Re b lr and Ira bi in equa- 
tions (75), (80), and (81), respectively, gives 

i re t + 2* y 

2^Ja t ^f(e ie )d9=log e ^ ( 86 ) 

1 r« ( -i-2x 

— Ref(e tB ) cos 6 dd— n cos 8, (87) 

1 re t -i-ar / 4 -TC , 2 \ 

-J a Rej{e ie ) sin 8 dd= sin 8, y 1 ^ (88) 

Hence, Re /(e fs ) must satisfy equations (86) to (88). 

Determination of Ref(e is ). — By matching the potentials 
4>c(s) and <f> t {6), a correspondence is established between 
points on the airfoil and the circle angles. The prescribed 
velocity distribution along the airfoil can therefore be written 
as a function of the circle angle q=q(B). “When absolute 
values of equation (67) are taken, 


Mg) \F'(e if )\ 

1 + Vl + iV) 


(89) 


where 


m 


Mg> 

l+Vi+a*® 


The values of Ref(e iS ), as given by equation (90), must 
satisfy equations (86) to (88). If these equations are not 
satisfied, the values of Ref{e iS ) must be adjusted until the 
equations are satisfied. One method for adjusting the 
function is to define 


Ref(e te )=Ref(e tB )—(jte b 0 —log e —(Re b y — ncos 8, ) cos 6— 

| ~Im bi— sin 8 t (n — ^ J sin 8 (91) 

where Re b 0 , Re b L , and Im b l are given by equations (83), 
(84), and (85), respectively. The function Ref(e ie ) will 
satisfy equations (86) to (88). If the corrections are large, 
this method may result in large changes in the velocity. 
In general, the method used for adjusting Ref(e iB ) will 
depend on the specific problem. 

Determination of Imf(e ie ) . — When ReJ{e u ) satisfying equa- 
tions (86) to (88) is obtained, Im f(e ie ) is computed by 
Poisson’s integral 

/ m /( g *9)=_ R e f( e if) C0i __ldT+( ai +ir) (92) 

where the constant term in the integral has been taken 
equal to ai+ir, so that 

* 

1 re,+2x 

7m6 0 =— Imf(e lS ) dd=cci-\-TT 

as required by equation (75). Hence j{e i0 ) is now known 
for points on the circle : 

J{e ie ) = Ref(e is )+ilmf(e is ) 


Determination ofg(e 1f ).— By use of these values of j(e te ) 
in equation (73), g(f) is determined for points on the circle: 

g(e is ) = [ l - e U»«-«>] v w ‘'' (93) 


for points on the circle. By replacing F'{e iB ) by the velocity 
v(6) (equation (69)) and substituting g(£) from equation (73) 
with equation (89) becomes 


2g(fl) 


2 1 7 sin 0+ ■^- L 
2x1 


l + Vl+ff'M) [9—2 cos (0«— fl)]*e a * /wM) 

or, with the equation solved for Ref(e m ), 


Rej{e i6 )—\og t 


2F sin 9+y 


[2 — 2 cos (8 t —8)] 2 K(8) 


( 90 ) 


Because of the changes made in Ref(e iB ) to satisfy equations 
-(86) to (88), g(e ie ) may no longer satisfy condition (65d). 
If g{e ts ) does not satisfy the inequality, g(£) can be adjusted 
by changing the second or higher harmonic terms in f(e ie ), 
or by some similar method. 

AIRFOIL COORDINATES 

The function gfy) has been obtained to satisfy all require- 
ments; hence, the airfoil coordinates in the 2 -plane are given 
by equation (64) on integrating around the unit circle. For 
convenience, let g(e ie ) be written as 

I g(e is )=g l (8)e i ^ 


056646— Bl- 
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From equation (64) the airfoil coordinates are then 


>-/ 


4fli 2 -(2Fsin 9+^J _ 




v ~i W, 


sin (jf 2 +0) dd (94) 


cos {g 2 + 6) d 6 (95) 


COMPUTATIONAL PROCEDURE FOR ISOLATED AIRFOIL 


An outline of the procedure for computing the airfoil 
follows: 

1. Obtain 4> c {s), r 0 , and <f> c (s n ) from equations (9) and (10). 
By use of these values in equations (71) and (72), obtain 
8 t and V, respectively. Then calculate and v(8) from 
equations (68) and (69), respectively. 

2. Plot 4>c(s) and 4> t {8). By comparing equal values of 
these potentials, obtain s as a function of 6, which permits 
writing the prescribed velocity q as a function of 8, q=q(8). 

3. Compute Ref(e ie ) by equation (90); this function must 
satisfy equations (86) to (88). If these equations are not 
satisfied, Ref(e ie ) is adjusted to satisfy them as suggested by 
equation (91) or by similar methods. 

4. Obtain Im j(e i0 ) from equation (92) and compute 
g(e ie ) by equation (93). These values of g(e ie ) should be 
checked in inequality (65d). If g(e ie ) does not satisfy the 
inequality, it is further modified as previously suggested. 

5. By use in equations (94) and (95) of the values of 
g(e iS ) that satisfy all requirements, the airfoil coordinates are 
obtained by integration. 


ILLUSTRATIVE EXAMPLES 


Isolated airfoil. — For this example, the prescribed velocity 
q=q(s) was obtained from the incompressible velocity dis- 
tribution about a symmetrical Joukowski profile, as com- 
puted by Lipman Bers at Syracuse University in the form 
of the ratio of actual velocity to free-stream velocity, by 
taking the dimensionless free-stream velocity to be 0.538. 
The resulting distribution is shown in figure 4 with the final 
velocity after adjusting Re jf(e w ). The computed profile 
(fig. 5) is slightly thicker than the Joukowski profile toward 
the nose and has some reflex camber. The peaks in the 
velocity distribution about the computed profile are lower 
than those that would occur in a compressible flow about the - 



Joukowski profile with the same free-stream velocity and 
angle of attack because: (1) the circulation was kept the 
same as for the incompressible flow, which resulted in the 
reflex camber; and (2) the thickening of the profile reduced 
the curvatures in the vicinity of the velocity peaks. 



o Joukowski Profile 
— ' — Computed Profile 

Figure 5.— Comparison of computed profile and Joukowski profile. 


Cascade. — For this example, it was decided to design a 
cascade of blades having cusped tails and a velocity distri- 
ibution along the blade like the distribution in the foregoing 
isolated-airfoil example. The free-stream conditions q u a u 
and a 2 were arbitrarily taken as 

gi=0.576 

oi=10° 

a 2 =0° 


which gives 


g 2 =0.564 

Adjustments to ReH(e te ) to satisfy equations (86) to (88) 
modified the prescribed velocity somewhat and the final 
velocity is shown in figure 6. Figure 7 is the computed 
cascade. 





q c *0.564 



Figure 4.— Comparison of prescribed velocity and final velocity distributions for airfoil. 


Figure 7. Two blades of computed cascade. 
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CONCLUDING REMARKS 

The magnitude of the dimensionless velocity along the 
blade cannot be entirely prescribed arbitrarily as a function 
of the arc length, but is subject to some restrictions in addi- 
tion to the conditions imposed on the regular functions H({) 
and /(f). The magnitude must be finite everywhere along 
the profile and, by the method given, the velocity can be 
zero in, at most, two places — the leading-edge and the 
trailing-edge stagnation points. By a limiting process, 
however, the method can be extended to provide for addi- 
tional stagnation points. The order of the zero of the dimen- 
sionless velocity at the trailing edge is determined by the 
included trailing-edge angle of the blade. Thus, for a cusp 
at the tail, the angle is zero and the dimensionless velocity 
need not be zero at the trailing edge. 


If a velocity distribution is selected to satisfy these con- 
ditions, but otherwise is arbitrary, the resulting profile may 
not be a physically real blade but may result in a blade with 
zero or negative thickness in some portions of the blade. 
The zero or negative thickness is caused by specifying too-low 
velocities along parts of the blade; a physically real blade 
can be obtained by increasing the prescribed velocity along 
the blade. 


National Advisory Committee for Aeronautics, 
Lewis Flight Propulsion Laboratory, 
Cleveland, Ohio, October 1, 1943. 


APPENDIX 


SYMBOLS 


The following symbols are used in this report: 

A,B, C i, Co, D complex constants 
«i,a 2 , . . . location of complex sources in f -plane 
C(f) function of f defined by equation (44) 

d spacing of cascade 

F(f) complex potential function (incompressible 

flow) 

/(f) regular function of f 

gr(f) regular function of f 

H( f) regular function of f 

Re H( f) function of f defined by equation (59) 


Im 

me) 

m 


m 

k 

n 

V 

Q(s) 

9 


q-,e ia * 

R 

Re 


V 
i’ IS) 

z^x+iy 


imaginary part 

function of B defined by equation (49) 
constant equal to — - ~ — . = - 

l + VT+i? 


constant equal to 


2g2 

1 + -/I+5? 


constant defined by equations (33) 
number determined by included trailing-edge 
angle of blade 
pressure 

auxiliary function of s 

magnitude of dimensionless velocity in 
compressible-flow plane (ratio of actual 
velocity to stagnation velocity of sound) 
dimensionless velocity upstream of cascade 
dimensionless velocity downstream of cascade 
region in f-plane defined by |f|> 1 
real part 

number defined by equation (23) 
arc length along blade 
free-stream velocity (incompressible flow) 
velocity on unit circle (incompressible flow) 
complex variable (compressible-flow plane) 


a 

r 

7 

5 

f=H-fi? 

6 
A 

p 

r 

0 

f 

Subscripts: 

c 

i 

n 

t 


angle of velocity in compressible flow (meas- 
ured from positive x-axis) 
circulation 

ratio of specific heats 
included trading-edge angle of blade 
complex variable (incompressible-flow plane) 
circle angle (incompressible-flow plane) 
auxiliary variable defined by equation (37) 
density 

variable of integration 
velocity potential 
stream function 


compressible flow 
incompressible flow 
leading edge 
trailing edge 


Prime indicates a derivative. 
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